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Abstract 

We study Fubini instantons of a self-gravitating scalar field. The Fubini instanton 
describes the decay of a vacuum state under tunneling instead of rolling in the presence 
of a tachyonic potential. The tunneling occurs from the maximum of the potential, which 
is a vacuum state, to any arbitrary state, belonging to the tunneling without any barrier. 
We consider two different types of the tachyonic potential. One has only a quartic term. 
The other has both the quartic and quadratic terms. We show that, there exist several 
kinds of new 0(4)-symmetric Fubini instanton solution, which are possible only if gravity 
is taken into account. One of them can be interpreted as the pair creation of an instanton 
in the de Sitter background. 
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1 Introduction 



The very first picture of an inflationary multiverse scenario was proposed in Ref. [T], in which, it 
would seem that the author wanted to suggest a universe without the cosmological singularity 
problem using an interesting feature of self-repro ducting or regenerating exponential expan- 
sion of the universe. A major development in this scenario was triggered by the discovery of 
the eternal inflationary scenario [2H5] and a paradigm for string theory landscape [6l[7]. The 
eternal inflation is related to the expanding false vacuum solution with a positive cosmological 
constant, which in turn means that the inflation is eternal into the future. If the theory has 
multiple minima then the false vacuum state decays into the true vacuum state, i.e. the phase 
transition proceeded from the nucleation of a vacuum bubble. In this scenario the universe 
is situated within some bubble called a pocket universe [1] having a certain value of the cos- 
mological constant and the whole universes are referred to as multiverse. The description of 
self-reproduction including tunneling process and random walk was combined into a scenario 
called recycling universe [S] . These scenarios seem to provide an escape from the question of the 
initial conditions of the universe, i.e. it seems to be eternal into the past. Unfortunately, infla- 
tionary spacetimes cannot be made complete in the past direction |9], even though the universe 
is eternal into the future. Anyway, this is not the central issue for our purpose and therefore we 
will not discuss this any further in the present paper. The string theory landscape is a setting 
that involves a huge number of different metastable and stable vacua pUjlTT]. originated from 
different choices of Calabi-Yau manifolds and generalized magnetic fluxes. The huge number 
of different vacua can be approximated by the potential of a scalar field. The important thing 
is the fact that, once the de Sitter vacuum can exist, the inflationary expansion is eternal into 
the future and has the possibility of self-reproduction. 

On the other hand, there are theories of gauged d = 4, N = 8 supergravity having de 
Sitter(dS) solution, in which all SUSYs are spontaneously broken. It is well known that the 
dS solution corresponds to a M/sting theory solution with a non-compact 7- or 6-dimensional 
internal space, in which a small value of the cosmological constant stems from the 4-form flux. 
The simplest representative of these kind of theories has a tachyonic potential with the dS 
maxima [T2] - [T^ . The potential in the vicinity of the maximum reduces to a form having a 
quadratic term, that is not metastable but unstable. However, according to some authors, the 
time for collapse giving rise to the tachyonic potential can be much greater than the age of the 
universe for anthropic reasoning. If the curvature radius of the potential in the vicinity of the 
maximum is greater than that used in the above theory, then that will be all together different 
story. The supergravity analogue of the tachyonic potential could be constructed also by using 
an exact supergravity solution representing the Dp-Dp system [T5] . 

From the above scenarios, the study of the possibility of the tunneling process for the 
potential with stable and metastable vacua, or even tachyonic behavior has acquired renewed 
interest. In the present paper, we will study the tunneling process under a simple tachyonic 
potential governed by a quartic term both without the quadratic term and with the term as a 
toy model. To obtain the general solution including the effect of the backreaction, we solve the 
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coupled equations for the gravity and the scalar field simultaneously. Although the model has 
a tachyonic potential, it might still be an useful example to show how the tunneling process 
occurs in various shapes of the potential provided by the above scenarios. 

A quantum particle can penetrate through a finite potential barrier via the so-called barrier 
penetration. This process can be described fully by in terms of the Euclidean solution obey- 
ing appropriate boundary conditions. There exist two kinds of Euclidean solutions describing 
quantum tunneling phenomena. One corresponds to an instanton solution representing a stable 
pseudoparticle configuration characterized by the existence of a nontrivial topological charge. 
It does not change even if we continuously deform the field, as long as the boundary conditions 
remain the same. The instanton solution corresponds to the minimum of the Euclidean action 
to pass from the initial to final state. The solution, in case of a double well potential, describes 
a general shift in the ground state energy of the classical vacuum due to the presence of an 
additional potential well, then lifting the so-called classical degeneracy. The other is a bounce 
solution representing an unstable nontopological configuration that corresponds to a saddle 
point rather than a minimum of the Euclidean action. The second derivative of the Euclidean 
action around the bounce has one negative eigenvalue which leads to the imaginary part of the 
energy. The existence of the negative eigenvalue implies that the vacuum state is unstable, i.e. 
the state decays into other states [16] . 

The Euclidean solutions can mediate phase transitions, too. The phase transition describes 
the sudden change of a physical system from one state to another. The transition are of two 
different types transition accompanied by temperature or without temperature. The competi- 
tion between the entropy and the energy terms in the thermodynamic potential cause thermal 
phase transitions in which dynamics is irrelevant. In the modern classification scheme, thermal 
phase transitions are divided into two broad categories either with a discontinuous jump in the 
first-order derivatives of the free energy or without it. A first-order phase transition is charac- 
terized by the discontinuity in the first derivative of the free energy and is associated with the 
existence of latent heat, whereas a nth-order phase transition is characterized by the continuity 
in the first derivative while there is a discontinuity in the nth-order derivative. A quantum 
phase transition describes a transition between different phases by quantum fluctuation, which 
occurs at zero temperature, unlike the case of a thermal phase transition which is governed by 
a thermal fluctuation [T7] . 

To simplify things, we consider an asymmetric double well potential to distinguish two 
different phase transitions without temperature. If the initial state is the metastable vacuum 
state and the tunneling occurs from that state to the other vacuum state, then the transition 
corresponds to a tunneling process [T9H23]. On the other hand, if the initial state is the local 
maxima of the potential and the fleld is rolling down to one vacuum state continuously rather 
than any discontinuous jump, then the transition corresponds to the rolling. However, one 
more channel exists as tunneling and that corresponds to the one without a barrier. In this 
kind of transition, the initial state on the top of a potential can tunnel to the other state 
rather than rolling down the potential [T8 | [25 | |27]. There are two different kind of transitions 
in this case. One is the tunneling without a barrier representing the tunneling from the local 
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maximum of the potential to the vacuum state [251 - I28] . Recently, an analytic study on this 
type of solution was performed in |2B] • The other is a tunneling without a barrier representing 
the tunneling from the maximum of the potential to any arbitrary state. This case corresponds 
to the Fubini instanton [18l|28], where the tachyonic potential is employed. Can we describe 
the rolling corresponding to the transition between the initial metastable vacuum state and the 
other final vacuum state? This may look similar to a superfiuid motion by the liquid helium. To 
establish the phase transition corresponding to the superfiuid motion is itself a very challenging 
problem, and hence we start by concentrating on the Fubini instantons in this work. 

The Fubini instanton [18] describes the decay of a vacuum state corresponding to a tun- 
neling process as the quantum phase transition instead of rolling down the tachyonic potential 
consisting of a quartic term only, in which the vacuum state is stable under small perturbations 
with finite energy. On the other hand, if one consider a tachyonic potential consisting of a 
quadratic term only, the point $ = is unstable. A small perturbation will cause it to roll 
down the hill of the potential. Originally, it was Fubini who introduced a fundamental scale of 
hadron phenomena by means of the dilatation noninvariant vacuum state in the framework of 
a scale invariant Lagrangian field theory [IB] . However, the solution is a one-parameter family 
of instanton solutions representing a tunneling without a barrier as an interpolating solution 
from the maximum of the potential to any arbitrary state. The instanton solution was studied 
in a conformally invariant model, i.e. a fixed background was used and the effect of the backre- 
action by instantons was neglected [25HHTj . This is a good approximation, when the variation 
of the potential during the transition is much smaller than the maximum of the potential. The 
instanton has gained much interest now-a-days in the context of anti-de Sitter(AdS)/conformal 
field theory correspondence [32|[33]. 

The paper is organized as follows: In Sec. II, we review the Fubini instanton in the absence of 
gravity. We present numerical solutions including the Euclidean energy density as an example 
and analyze the structure of the solution in the theory with a potential having only the quartic 
self-interaction term. We stress the fact that there is no such solutions with the potential 
containing both the quartic and the quadratic terms. In Sec. Ill, we show that the instanton 
solutions exist in the curved space. We perform a numerical study to solve the coupled equations 
for the gravity and the scalar field simultaneously. We show that there exist numerical solutions 
without oscillation in the initial AdS space in the potential without any quadratic term. We 
also show that there exist numerical solutions in the potential both with the quartic and the 
quadratic terms irrespective of the value of the cosmological constant, which is possible only 
when the gravity is switched on. In order to estimate the decay rate of the background state, 
we compute the action difference between that of the solution and the background obtained 
by numerical means. In Sec. IV we present an oscillating numerical solutions in the potential 
without any quadratic term regardless of the value of the cosmological constant. One type of 
these solutions can be interpreted as the pair creation of a Fubini instanton in the de Sitter 
background. We analyze the behavior of the solutions using the phase diagram method. To 
observe the dynamics of the solutions, we briefly sketch the causal structures of the solutions 
in the Lorentzian spacetime. Finally in Sec. V, we summarize and discuss our results. 
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2 Fubini instanton in the absence of gravity 



One can consider the following action in the absence of gravity 



S = X 
'm 



where g = detrjn^, and the tachyonic potential has a quartic self-interaction term and also a 
quadratic term as follows: 

[/($) = + + f/„, (2) 

where > and A > 0. The plots of potentials (a) without the quadratic term and (b) 
with the quadratic term are shown in Fig. [1] The potential has a metastable vacuum state at 
$ = and no other stationary state in Fig. [U^a), while Fig. HJ^b) illustrates that the potential 
has a local minimum at $ = and two maxima . In both the cases, the potential is not 

bounded from below. 

Before going to the tunneling problem in four dimensions, we briefly describe the problem 
in one dimension. One can consider the simplest quantum tunneling problem in one dimension. 
Quantum field theory in one dimension is nothing but ordinary quantum mechanics. In case 
of = 0, the amplitude for transmission obeys the WKB formula in the semiclassical ap- 
proximation, in which $-1- = ±(^)^/'* are the classical turning points. On the other hand, the 
double-hump potential with > and Uo = can be considered as an inverted double-well 
potential for a bounce solution representing the tunneling from $ = to $-1- = ±m-^2/A. The 
solution is given by $s(t) = ±m^2/Asech [m(r — Tq)], where To is an integration constant. The 
bounce solutions can be easily understood to Euclidean space. The particle can only reach 
the point $ = Oatr = ±oo and it bounces off $-1- = ±m^y2/X at r = with a vanishing 
velocity [37] . 

We now turn to the tunneling problem in four dimensions. It is an well-known fact that 
the massless theory has an instanton [18]. Actually, the instanton corresponds to the bounce 
solution representing the decay of the background vacuum state. The equation of motion with 
0(4) symmetry, obtained by varying the Euclidean action, is then; 

c/2$ ^ 3 d$ _ d{-U) 
dr]"^ 7] dr] d^ ' 

where the second term in the left-hand side plays the role of a damping term. The boundary 
conditions are 

— =0 and <|.|,,=^ = 0. (4) 
dr] 77=0 

The particle in the classical mechanics problem starts at $ = with zero velocity in the 
inverted potential, and stops at $|r,=oo = without any oscillation. 
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Figure 1: Potentials for the case of (a) Fubiiii instantons and (b) generalized Fubini instantons. 
For the potential with = 0, the analytic solution of the Fubini instanton has the form 

where r] is the radial length in the Euclidean space, b is any arbitrary length scale which 
characterizes the size of the instanton and is related to the initial value In addition, the 



value of the scalar field of the center of the solution depends on b as $(0) = This solution 

was widely used in the formulation of Yang-Mills instantons [38] and the related perturbation 
theory [39] . 

The characteristic behavior of the analytic solution is ploted in Fig. [2] in terms of the value 
of the parameter b. We take A = 1 for all the cases. The solid line denotes the solution with 
b = 1, the dashed line with b = 3, and the dotted line with b = 5. The corresponding Euclidean 
action is given by 



where the action does not depend on the parameter b due to the consequence of the conformal 
invariance of the potential and we take that value to be Uo = 0. The action has the same value 
irrespective of the starting point $0. In other words, the tunneling from the maximum of the 
potential to any arbitrary state always happens with same probability. 

The numerical solutions for $ and $' and the Euclidean energy including the density varia- 
tion with rj are as shown in Fig. [31 Figure [31^a) illustrates the numerical solution for $, in which 
the initial value set as $0 = — 1 and the solution asymptotically approaches the value $ = 0. 
Figure [3l^b) illustrates $' with respect to rj. There is a peak of $' near t] = 2.31. Figure |3t^c) 
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Figure 2: The analytic solution of the Fubini instanton in absence of gravity. 

depicts the volume energy density, when the density has got the form ^ = [|$'^ + U]. The lower 
right box in the same figure shows the magnification of the small region clearly representing the 
existence of a smooth hill. The smooth peak of the volume energy density exists at 77 = 5.17. 
There is a disagreement between the location of the peak for the energy density and that for 
It clearly reveals the fact that the position with the maximum value for $' is still not the 
same as the maximum of the energy density due to non-trivial contribution coming from the 
potential, U = —j^"^. Figure [3]^d) shows the Euclidean energy for each slice of constant 77. 
The Euclidean energy signifies the value of energy after the full integration of variables except 
for r] in the present work, = 2Ti'^7f^. There are one minimum and one maximum point for 
E(^. The location of the minimum of E^ is around rj = 2.31, whereas that of the maximum 
is near rj = 6.93. Ironically, the location of the minimum of E^ coincides with that of the 
maximum of These solutions can be considered as a ball consisting of only a thick wall 
except for one point at the center of the solution with a lower arbitrary state than the outer 
vacuum state unlike a vacuum bubble that consists of an inside part with a lower vacuum state 
and a wall. 

For a theory with > 0, the conformal invariance is broken and a regular solution with a 
finite action is forbidden by scaling argument, i.e. there exists no instanton solution with 0(4) 
symmetry. In other worlds, the particle can not have enough energy to reach the hill overcome 
the barrier near $ = since the damping term has got a large value due to a large value of $' 
near the initial point |40j . 
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Figure 3: (a) The numerical solution for $ in the case of = 0, (b) the variation of $' with respect to t], 
(c) the energy density ^ and (d) the Euclidean energy i?^ evaluated at constant 77. 



3 Fubini instantons of a self- gravitating scalar field 



Let us consider the following action: 



S 



M 



-h(fx- 



(7) 



dM 



where g = detf^^jy, n = SirG, R denotes the scalar curvature of the spacetime Ai, h is the 
induced boundary metric, K and Ko are the traces of the extrinsic curvatures of dAi for 
the metric g^u and rj^^, respectively. The second term on the right-hand side is the boundary 
term ^1] . It is necessary to have a well-posed variational problem including the Einstein-Hilbert 
term. Here we adopt the notations and sign conventions of Misner, Thorne and Wheeler |42j . 

We study the creation process of Fubini instantons in curved spacetime. In the first place, 
we consider the massless case, and then we will also consider generalized Fubini instantons, the 
so-called massive case (see the form of the potential in Eq. ([2])). The cosmo logical constant is 
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given by A = kUq, such that background space will be dS, flat or AdS depending on the signs 

In order to solve the coupled equations, we assume an 0(4) symmetry for the geometry and 
the scalar fleld similar to Ref. [2T] 



ds^ = drf + p2(^) [dx^ + sin^ X {d6^ + sin^ Odcj?)] . (8) 

And then, $ and p depends only on r/, and the Euclidean equation can be written respectively 
as follows: 

$" + = ^ and p" = + U) , (9) 

p rf$ 3 

and the Hamiltonian constraint is then given by 

p'2_i_1P^[l^'2_u\ (10) 



3 V2 

In order to yield a meaningful solution, the constraint requires a delicate balance among all the 
different terms. Otherwise the solution can yield qualitatively incorrect behavior |13]. 

To solve the Eqs. ([9]), we have to impose suitable boundary conditions. When the gravity 



is switched off, boundary conditions for the Fubini instanton are , 
in Ref. [18]. While gravity is taken into account, the initial conditions are as follows: 



and $|,,=oo = as 



I n 



= 1, $1^=0 = and — 

r;=0 drj 



= 0, (11) 



where the flrst condition is to obtain a geodesically complete spacetime. The second condition 
is nothing but Eq. f llOp . The fourth condition is the regularity condition as can be seen from the 
flrst equation in Eq. (|9]). The third condition for the initial value of $ remains undetermined. 
One should flnd the initial value of $ using the undershoot-overshoot procedure [20l[27] . 

In order to solve the Euclidean fleld Eqs. ([9]) and ( ITO|) numerically, we rewrite the equations 
in terms of dimensionless variables as in Ref. [27] . In the present work, we employ the shooting 
method using the adaptive step size Runge-Kutta as an initial value problem similar to the 
treatment in Ref. [H]. For this procedure we choose the initial values of $(?7initiai)) '^''(^initiai), 
p(?7initiai), and p'(r/initiai) at 77 = ?7initiai as follows: 



(12) 



(^initial) 


~ - -|.,(<|2 - 1 


(^initial) 


1<|.,(<|2 _ 1) + . 


P(^initial) 


~ e + ■ ■ ■ , 


P'(r7initial) 


~ 1 + ■ ■ ■ , 



where ?7initiai = + e for e ^ 1. The minus sign in front of the second formula is due to the 
negative value of the $" determined by the sign of dU /d^ at = 0. However, the initial value of 
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Table 1 : The dimensionless variables and color of plot used and the actions obtained in Fig. 21 

$' is taking to be positive in the present work. Once we find the initial value $o, the remaining 
conditions can be exactly determined from Eqs. ( IT2l) . The parameter k is the ratio between the 
gravitational constant or Planck mass and the mass scale in the theory, k = = ^py, and 

the parameter kUo is related to the rescaled cosmological constant A./m?. 

If the size of the solution is finite we should impose same additional conditions implicitly 
when we consider the pair creation of the instantons. To avoid a singular solution aX fj = fjma.x 
for the Euclidean field equations and to demand a Z2 symmetry, the conditions d^/dfj — and 
p ^ as fj ^/max are needed in the next section. In this work, we require that the value of 
d^/dfj goes to a value smaller than 10~^ as — )■ r^max, as the exact value of r^max is not known. 

To find the probability of the instanton solution, we only consider the Euclidean action for 
the bulk part in Eq. ([7]) to get. 



I Vg^d'xE + + u\= 2n' [ p^dv[-U] , (13) 

Jm I Zk I J J 



where -Re = 6[l/p^ — p'"^ / — p" / p\- We used Eqs. ([9]) and ( fTOj) to arrive at this. The volume 
energy density has the form: ^ = — f/, which has a different sign compared to the sign of 
the density used in Ref. [2^. The Euclidean energy signifies the energy value after the full 
integration of variables except for 77 in the present case as = 2-7r^p'^^. 

In the beginning, we obtain the numerical solution for = 0. And then we obtain the 
numerical solution for > 0. We call the space dS when the initial vacuum state has a 
positive cosmological constant, Uo > 0, flat when Uo = and AdS when Uo < 0. 

The rate of decay of a metastable state can be evaluated in terms of the classical con- 
figuration and represented as Ae~^ in this approximation, in which the leading semiclassical 
exponent B = S'^^ — S^^ is the difference between the Euclidean action corresponding to the 
classical solution S'^^ and the background action S^^. The prefactor A is evaluated from the 
Gaussian integral over fluctuations around the background classical solution j45] . 

The numerical solutions for instantons with = in the AdS space are shown in Fig. |H We 
take K = 0.1. The solutions presented in this paper are possible in the initial AdS background 
only for the case with = 0. We are not quite sure whether the solution without any 
oscillation is still possible for the initial fiat and dS background or not. Figure 4(a) illustrates 
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Figure 4: (color online). The numerical solutions of Fubini instantons with 



in the AdS space. 



the solution for $, in which the right box in the same figure shows the magnification of a 
small region representing the initial values of $ and the behavior of the curves. The curves 
move upwards with increasing value of [/q, then overlap near r/ = 1, and more downwards with 
increasing value of Vo- Figure 4(b) shows the solutions of p. The curves move downwards with 
increasing Vo- The shape of the numerical solution p can be easily understood if one thinks 

of the shape of the solution in a fixed AdS space as p = y^sinh -y/f ^7- Table 1 shows the 
dimensionless variables and the color of plot used, and also the actions obtained from Fig. HI 
From the numerical data, one can easily see that the magnitude of $o approaches the vacuum 
state $ = as f/o approaches a vanishing value. The vanishing of Vo means the background 
geometry which serves as the initial vacuum state is fiat. The action difference B between the 
action of the solution S""^ and that of the background S^^^ has positive values. We carry out 
the action integral in the range ^ ^ 25 numerically as the action difference B diverges 
to infinity if we perform the integration for an infinite i] value. This divergence is due to the 
fact that the size of the solution including the outside part in the evolution parameter space 
decrease compared to the size of the initial background similar to what happens for the case 
of the nucleation of a vacuum bubble. In the analytic computation, the outside part and the 
background are simply canceled at the same radius. In the present numerical work, it is difficult 
to decide the exact size of the solution. Thus we straightforwardly compute the action difference 
and then the difference B has got an approximate behavior 5(sinh^ r/) = 3 sinh^ 77 cosh 77 which 
cause the divergence at infinity. If this minor error is cured, the action difference has a finite 
value. 

The numerical solutions for the case with rn? > in curved space with an arbitrary cos- 
mological constant are shown in Fig. |5l We take k = 0.3 for this case. The figures represent 
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Figure 5: (color online). The numerical solution of $, the derivative of $ with respect to rj, p, and the 
derivative of p with respect to rj for the generalized Fubini instantons with m? > 0. 

the vary fact that the solution is only possible in curved spacetime irrespective of the value of 
the cosmo logical constant. Figure 5(a) illustrates the solution of $. The upper right box in 
the same figure shows the magnification of the small region representing behavior of the curves 
which move to the left with an increase in Uo- Figure |5l^b) shows $' with respect to rj. The upper 
right box in the same figure shows the magnification of the small region representing behavior 
of the curves moving below with increasing value of f/g. Figure EI^c) illustrates the solutions of 
p. The curves move downwards with increasing value of Uo- In the fixed fiat space, p = f]- In 

the dS space, p = y^sin \J^V- Figure |5t^d) depicts the variation of p with respect to rj. The 

curves move below with increase in Uo- The horizontal line with Uo = indicate a fiat space 
with p' = 1- Table |2] shows the dimensionless variables and color of the plot used among with 
the action obtained from Fig. |5l From the numerical data, one can infer that the magnitude 
of $o decreases as Uo increases. We carry out the action integral in the range ^ fj ^ 30.58 
numerically. In the dS space, the solution and the background have their own periods for fj, 

which we as the integration limit. For the background dS space, we take fj = tt./t^. The 
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Table 2: The dimensionless variables and color of plot used and the actions obtained in Fig. [5l 

action for S'^'^ and S^^ are positive or zero as long as Uo ^ 0. The background action is zero 
for Uo = 0. In this work, we do not check for the special case S'^^ = 0. Simply, the action 
has a negative value for the dS space. It is related to the fact that the Euclidean action for 
Einstein gravity is not bounded from below, and this is known as the conformal factor problem 
in Euclidean quantum gravity [16]. It was argued in ^\ that the conformal divergence due 
to the unboundedness of the action might get cancelled with a similar term of opposite sign 
caused by the measure of the path integral. However, the difference between the action of the 
solution and that of the background remains positive-valued. 

4 Oscillating Fubini instantons 

The oscillating instanton and the bounce solutions with an 0(4) symmetry between the dS-dS 
vacuum states was first studied in Ref. [IS], in which the authors found the solutions in a fixed 
background geometry and showed how does the maximum allowed number Umax depend on 
the parameters of the theory, where n denotes the crossing number of the potential barrier by 
the oscillating solutions. The oscillating bounce solutions in the presence of gravity was also 
studied in Ref. [19], in which the authors analyzed the negative modes and the fluctuations 
around the oscillating solutions. The oscillating instanton solutions under a symmetric double- 
well potential in the curved space with an arbitrary vacuum energy was also investigated in 
detail in [27], where a numerical solution is possible as long as the local maximum value of 
the potential remains positive. The solutions have a thick wall and can be interpreted as a 
mechanism for the nucleation of the thick wall for topological inflation [50]. Similarly, the 
process for the tunneling without a barrier in curved space, was studied in Ref. [SSH^B]. The 
existence of numerical solutions was shown in Ref. [27], in which the case representing the 
tunneling from flat to AdS space shows an oscillating behavior. The solution oscillates around 
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Table 3: The dimensionless variables and color of plot used and the actions obtained in Fig. [Bj 

$ = in the inverted potential and the oscillating behavior die away unlike the case under a 
harmonic potential. 

In the present paper, the oscillation means that the field in the solutions oscillates around 
the minimum of the inverted potential and die away asymptotically to the minimum $ = 
for the case with = 0. Thus the resulting geometry of the initial state has wrinkles due to 
the variation of the volume energy density and the instanton simultaneously. The behavior of 
the solutions representing the resulting geometry with wrinkles is quite different from those in 
Ref. [27]. 

Figure [6] shows the numerical solutions representing an oscillatory behavior in (A) the initial 
dS background, (B) the initial flat background and (C) the initial AdS background. We take 
R = 0.30, Uo = 0.0001 (for the dS case), Uo = (for the flat case), and Uo = -0.0001 (for 
the AdS case), respectively. Figures |6t^ a), (c) and (e) illustrate the numerical solutions of $. 
The lower right box in those figure shows the magnification of a small region representing the 
behavior of the solution around $ = 0. The peak corresponds to the first turning point of the 
particle similar to a classical mechanics problem in the presence of an inverted potential. For 
the case with $o = —7 (yellow line) the first turning point reaches furthermost point away from 
$ = among all the other cases, as one can easily see from the figure. The curves oscillate 
around $ = and eventually stop at $ = in the fiat and AdS space. We also adopt a similar 
behavior representing the stopping of the curve at $ = also in dS space. Figures El^b), (d) 
and (f) illustrate the numerical solutions for p. The upper left box in those figure shows the 
magnification of an initial small region representing behavior of the curves which move below 
with the decrease in $o. Table [3] shows the initial values of $, colors of the plot we used and 
the actions for the AdS and fiat background which are obtained from Fig. |6l In the fiat case, 
the background action is zero as f/o = and therefore S'^^ is equal to B. In the dS case, S'^^ is 
rapidly decreasing compared to S^^ which in turn means that the difference between the action 
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Figure 6: (color online). The numerical solutions representing oscillatory solutions. 
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of the solution and that of the background has negative value at early times but it will become 
positive at late time. In the present case, we cut all the data at a certain point which is f/ = 200 
and the evolution parameter of the dS background does not even fully exceed the given range. 
Thus it is hard to express the full action in a dS background. However we can expect that the 
action in the dS background is slightly smaller than that in the flat case as can be inferred from 
the result of Table [2l 

We now analyze the behavior of the solutions using a phase diagram method. After plugging 
the value of — from Eq. (fTOj) into Eq. ([9]) and using $" = the equation becomes 



First, we consider the situation where the kinetic energy is small compared to the potential 
energy such that ^ ?7o <C 1 and the term is smaller than other terms. In other 
words, the last term is the most dominant among other terms in the numerator. Then the 
equation reduces to the form 



The above relation shows that the first stage of the curve has got such kind of form. Second, 
we consider the situation where d^'/d^ = 0, i.e. with vanishing acceleration and then we 
impose all the above mentioned conditions. It will then describe the special points in the phase 
diagram. Then the equation reduces to the form 



Third, we consider the situation where d^'/d^ = —c, i.e. a negative constant. We impose 
all the above mentioned conditions among with $^ 2c/a/3kA. Thus, we obtain the above 
equation again. This relation implies that the special points with a vanishing acceleration and 
some of the region with a negative constant acceleration in the phase diagram have got a linear 
function type behavior in the phase diagram as shown in Figs. [7|(a), (c) and (e). 

Figure [7] illustrates the behavior of the solutions in the plane. Each trajectory rep- 
resents the behavior of the solution in the phase diagram. The trajectories begin with zero 
velocity as $' = shown in Figs. El^a), (c) and (e). The velocity increases rapidly to the max- 
imum and then decreases linearly up to the turning point. Figures [7](b), (d) and (f) show the 
magnification of the small region representing the behavior of the solution around $' = and 
$ = 0. 

Basically, the Fubini solution has an asymptotic condition to be satisfied. We expect that 
there exist an oscillating solutions although the dS background has got a finite size in the Eu- 
clidean signature. However, if we consider an analytic continuation not of the angle parameter 
X but of the Euclidean evolution parameter rj = it, then the meaning becomes clearer. When 
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Figure 7: (color online). The behavior of the solutions represented in the phase diagram. 
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there is an 'even' symmetry for the oscillating instantons, we can see the half-way point rjo as 
p(r7o) = $(?7o) = 0. Then, we can paste the Lorentzian manifold t = at the f] = rjQ surface. 
This is possible only for the case p(?7o) = ^('7o) = 0, because of the Cauchy-Riemann theorem 
of complex analysis; otherwise, the Lorentzian manifold should be complex valued functions 
(for exceptional cases, we might be able to consider complex valued instantons, the so-called 
fuzzy instantons [5TII52] ). In this procedure, an event shows a spontaneous creation of the 
universe from 'nothing' [2], in which nothing means a state without the concept of classical 
spacetime [53] • We already know that there is such a solution when the scalar field is ex- 
actly on top of the local maximum. However, now we observe a creation from nothing with 
highly non-trivial field dynamics. This is worthwhile to be highlighted and we postpone further 
analysis for the future work. 

4.1 Pair creation of Fubini instantons 

We now shift our attention to the pair creation of Fubini instantons. The pair creation of the 
instantons in this work is quite different from the ordinary quantum process of pair creation 
of particles. We take the initial background as the dS space, i.e. Uo > 0. The Euchdean dS 
space has a compact geometry. Thus, the geometry has two poles. If one object is created 
on the north pole and the other on the south pole, we can interpret that process as the pair 
creation of objects. As an example of the process, the two-crossing solution between the sides 
of the potential barrier in the double-well potential was considered as a type of a double-bounce 
solution or an anti-double-bounce solution [51], in which the authors interpreted the double- 
bounce solution as the spontaneous pair-creation of true vacuum bubbles, one at each pole in 
the dS space. We adopt a similar interpretation for our pair creation process. Figure [H] shows 
the numerical solutions representing pair creation of the Fubini instanton. We take k = 0.50 
and Uo = 0.03. Thus the dS region in the $-space spans the region —0.589 < $ < 0.589. We 
consider four cases with different initial positions of $. Figure |8]^a) illustrates the numerical 
solution of $. Figure MJo) depicts the numerical solution of p. Figures [S](c) and (d) illustrate 
the behavior of the solutions in the plane. Each trajectory represents the behavior of 
the solution in the phase diagrams. The blue and red lines indicate that the interior part of 
two instantons has got the same state as whereas the green and black lines indicate that 
the interior part has got the opposite state of $. Figure [Ht^d) illustrates the magnification of 
a small region representing the behavior of the solution around $' = and $ = 0. Figure 
illustrates the volume energy density, where the density has got a form ^ = —U. The box 
shows the magnification of a small region representing behavior of curves. The densities in each 
of the case have got positive values near the initial starting point $0 far away from the point 
$ = 0, because the densities have the form ^ = —U and Uo > 0. The solutions oscillate in the 
dS region as found the present work. The density always negative values for the case of the 
black line. Figure IHI^f) illustrates the Euclidean energy = 27r^p^^ for each slice of constant r] 
values. The negative energy parts in each of the case signifies a rolling state in the dS region. 
Table m shows the initial value of $, colors of the plot used, and the actions obtained from Fig. 
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Figure 8: (color online). The numerical solutions representing the pair creation of the Fubini instanton. 

El 

The trajectories with the blue and red color go back to the same position of $ in the 
presence of the inverted potential, whereas the trajectories with the black and green color go 
back to the opposite position of $. This can be interpreted as a pair creation with each one 
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Table 4: The dimcnsionless variables and the color of plot used, and the actions obtained in Fig. [Hj 

having the same state for the cases with the blue and red color, and each one having opposite 
state for the cases with the blue and red color, respectively. The solutions they can have 
more interesting interpretation as follows: any arbitrary state can tunnel to the same arbitrary 
state with an (9(4)-symmetry in curved spacetime, even though no vacuum state, as exists an 
instanton solution. The solutions are possible as long as the maximum of the potential remains 
positive. 

4.2 Causal structures 

In this section, we briefly outline the causal structure of the solutions in the Lorentzian signa- 
ture. Due to the pressure difference, the nucleated AdS region will expand over the background 
and hence the boundary of the nucleated AdS region will be time-like. 

Figure |9] shows the schematic diagrams representing the causal structures of the Fubini 
instantons and the related solutions. The x = 7r/2 surface can be analytically continued to the 
surface t = in the Lorentzian signature. The lower vacuum region in the instanton (green 
colored region) will be unstable during the Lorentzian time evolution (orange colored region). 
Due to the instability of the Fubini type potential, the whole causal structure may depend on 
the shape of the potential or the vacuum structure i.e. whether the left or the right side of the 
potential has true vacua or not. Therefore, the followings are meaningful only as reasonable 
estimations for general behavior and these may be different for special examples. 

Figure M^a) illustrates the instanton solution in an AdS background. It will form time-like 
r = and r = oo boundaries in the Lorentzian signature. However, the AdS region is always 
unstable to form a kind of singularity. Figure Mjo) illustrates the instanton solution in the 
dS background. The dS region has a cosmological horizon and will this form a future infinity. 
The AdS region (orange colored region) will expand over the dS region due to the pressure 
difference. Figure M^c) is the pair creation by the oscillating instanton solutions. Therefore, 
in the instanton part, the dS region around the p = Pmax is surrounded by the AdS (green 
colored region) part. In the Lorentzian signature, we can interpret these two AdS parts as 
being nucleated in a dS background. In Fig. M^c), we infer that, there still remains a dS region 
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(a) (b) (c) 

Figure 9: (color online). The schematic diagrams representing the causal structure of the Fubini instantons 
and the related solutions, (a) Tunneling in an AdS background, (b) oscillating instanton solution in a dS 
background, (c) Pair creation by oscillating instantons in dS background. 

and a future infinity; but it does depend on the choice of the potential parameters and if the 
instanton expands relatively quicker, then the AdS regions might occupy all of the dS region 
although this is not what happens in our case. 

5 Summary and Discussions 

In this paper we have studied Fubini instantons of a self-gravitating scalar field representing 
the tunneling without a barrier. There are two kinds of Euclidean solutions representing the 
tunneling without any barrier. One of them is the tunneling from the local maximum of the 
potential to the vacuum state. The other one is the tunneling from the maximum to any 
arbitrary state. The latter corresponds to the Fubini instanton solution. We have shown that 
there exist several new kinds of Fubini instanton solutions of a self-gravitating scalar field 
found as numerical solutions, which are possible only if gravity is taken into account. We also 
computed the action difference B, in each case, between the Euclidean action corresponding to 
a classical solution S'^'^ and the background action S^^ for the rate of decay. 

In Sec. II, we reviewed the Fubini instanton in the absence of gravity from the viewpoint of 
a tunneling problem. We have presented a numerical solution including the Euclidean energy 
density for example. We analyzed the structure of the solution in a theory with the potential 
having only a quartic self-interaction term. These solutions can be considered as a ball consisting 
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of only a thick wall except for the one point at the center of the solution with a lower arbitrary 
state than the outer vacuum state unlike a vacuum bubble which consists of an inner part with 
a lower vacuum state and a wall. 

In Sec. Ill, we have studied the instanton solutions in curved space. We performed careful 
numerical study to solve the coupled equations for the gravity and the scalar field simultane- 
ously. We have shown that there exist numerical solutions without any oscillation in the initial 
AdS space for the potential without the quadratic term. We have also shown that there exist 
numerical solutions for the potential with both a quartic and a quadratic term irrespective 
of the value of the cosmo logical constant. For this particular case, there is no solution with 
an 0(4) symmetry when gravity is switched off. In order to estimate the decay rate of the 
background state, we calculated the action difference between the action of the solution and 
that of the background obtained using numerical means. 

In Sec. IV, we have obtained oscillating Fubini instantons as new types of solutions. We have 
shown that there exist oscillating numerical solutions for the potential without any quadratic 
term regardless of the value of the cosmological constant. We have analyzed the behavior of 
the solutions using the phase diagram method. For the case of the dS space, there are some 
difficulties to interpret the solution. However, the oscillation dies away asymptotically in both 
the flat and the AdS space. Thus, we extend the scope of the oscillation to the dS space even 
including an Euclidean dS space having a flnite geometry. 

We have obtained numerical solutions representing the pair creation of Fubini instanton. 
These solutions can be interpreted as the pair creation with each one having the same state 
and with each one having the opposite state, respectively. The solutions can lead to more 
interesting interpretation as follows: any arbitrary state can tunnel into another arbitrary 
state with an 0(4)-symmetry in the curved spacetime, although no vacuum state exists as the 
instanton solution. The solutions are possible as long as the maximum of the potential remains 
positive. 

Finally, we have analyzed the schematic diagrams representing the causal structures of the 
Fubini instantons and the related solutions in the Lorentzian signature. Due to the pressure 
difference, the nucleated AdS region will expand over the background. For the special case 
representing the pair creation, the dS region around the p = Pmax is surrounded by an AdS 
part. 

In summary, we illustrate the following flnding in our new contribution regarding this issue: 

1. In the absence of gravity, a — </)^-type potential has inflnitely many instanton solutions 
whereas a —0^ + 0^-type potential has no instanton solution. However, the inclusion of 
the gravity changes all the situation abruptly: for the former case, the solution space get 
reduced to a unique solution and for the latter case, there exists solutions. 

2. We also conflrm that — 0^-type potentials have oscillating instanton solutions as well as 
the pair creation by the oscillating solutions. 

Therefore, the Fubini instanton is one of the few examples that shows the effect of gravity 
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bringing drastic changes to the tunnehng process. There can be more apphcations of the 
oscillating instantons and we confirm that the Fubini-type potentials also contribute largely 
towards these processes. We postpone any possible application of such oscillating solutions for 
our future work. 
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